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We investigate in detail the question whether a non-vanishing cosmological constant is required by 
present-day cosmic microwave background and large scale structure data when general isocurvature 
initial conditions are allowed for. We also discuss differences between the usual Bayesian and the 
frequentist approaches in data analysis. We show that the COBE-normalized matter power spectrum 
is dominated by the adiabatic mode and therefore breaks the degeneracy between initial conditions 
which is present in the cosmic microwave background anisotropies. We hnd that in a flat universe 
the Bayesian analysis requires Da 0 to more than Scr, while in the frequentist approach Da = 0 is 
still within Scr for a value oi h < 0.48. Both conclusions hold regardless of initial conditions. 

PACS numbers: PACS: 98.70Vc, 98.80Hw, 98.80Cq 


I. INTRODUCTION 


Ever since the beginning of modern cosmology, one of 
the most enigmatic ingredients has been the cosmological 
constant. Einstein introduced it to find static cosmolog¬ 
ical solutions (which are, however, unstable) 0]. Later, 
when the expansion of the Universe had been established, 
he called it his “greatest blunder”. 

In relativistic quantum field theory, for symmetry rea¬ 
sons the vacuum energy momentum tensor is of the form 
for some constant energy density e. The quantity 
A = SttGc can be interpreted as a cosmological constant. 
Typical values of e expected from particle physics come, 
e.g., from the super-symmetry breaking scale which is 
expected to be of the order of e>l TeV"^ leading to 
A>1.7 X 10“^® GeV^, and corresponding to Da^IO^®. 
Here we have introduced the density parameter Da = 

e/Pcrit = A/(87rG/9crit), where pcrit = 8.1 x 10“'^'^/i^ GeV^ 
is the critical density and the fudge factor h is defined 
by Hq = 100 h km s“^ Mpc“^, it lies in the interval 
0.5 < h < 0.8. Hq is the Hubble parameter today. 

Such a result is clearly in contradiction with kinemat- 
ical observations of the expansion of the universe, which 
tell us that the value of Dtot, the density parameter for 
the total matter-energy content of the universe, is of the 
order of unity, ©(Dtot) ~ 1- For a long time, this appar¬ 
ent contradiction has been accepted by most cosmologists 
and particle physicists, with the hind thought that there 
must be some deep, not yet understood reason that vac¬ 
uum energy — which is not felt by gauge-interactions — 
does not affect the gravitational field either, and hence 
we measure effectively A = 0. 

This slightly unsatisfactory situation became really 
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disturbing a couple of years ago, as two groups, which 
had measured luminosity distances to type la supernovae, 
independently announced that the expansion of the uni¬ 
verse is accelerated in the way expected in a universe 
dominated by a cosmological constant The ob¬ 

tained values are of the order 0(D„i) ~ 0 (Da) ^ 1 and 
cannot be explained by any sensible high energy physics 
model. Tracking scalar fields or quintessence Q ||] and 
other similar ideas Q have been introduced in order 
to mitigate the smallness problem — i.e., the fact that 
e ^ 10“"^® GeV'*. However, none of those is completely 
successful and really convincing at the moment. 

After the supernovae la results, cosmologists have 
found many other data-sets which also require a non¬ 
vanishing cosmological constant. The most prominent 
fact is that anisotropies in the cosmic microwave back¬ 
ground (GMB) indicate a flat universe, Dtot = Di„+Da = 
1, while measurements of clustering of matter, e.g., the 
galaxy power spectrum, require T = hflm ~ 0.2. But also 
CMB data alone, with some reasonable lower limit on the 
Hubble parameter, like h > 0.6, have been reported to 
require Da > 0 at high significance (see, e.g., Refs. |^, ||] 
and others). 

This cosmological constant problem is probably the 
greatest enigma in present cosmology. The supernova re¬ 
sults are therefore under detailed scrutiny. For example 
Ref. is not convinced that the data can only be un¬ 
derstood by a non-vanishing cosmological constant. Gos- 
mological observations are usually very sensitive to sys¬ 
tematic errors which are often very difhcult to discover. 
Therefore, in cosmology an observational result is usu¬ 
ally accepted by the scientific community only if several 
independent data-sets lead to the same conclusion. But 
this seems exactly to be the case for the cosmological 
constant. 

This situation prompted us to investigate in detail 
whether present structure formation data does require 
a cosmological constant. One may ask whether enlarging 
the space of models for structure formation does miti- 
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gate the cosmological constant problem. There are sev¬ 
eral ways to enlarge the model space, e.g. one may al¬ 
low for features in the primordial power spectrum, like a 
kink 0 - In the present paper we study the cosmologi¬ 
cal constant problem in relation to the initial conditions 
for the cosmological perturbations. In a first step we re¬ 
discuss the usual results obtained assuming purely adi¬ 
abatic models and we investigate to which extent CMB 
alone or CMB and large scale structure (LSS) require 
Ha 0 in a flat universe. We shall hrst present the 
usual Bayesian analysis, but we also discuss the results 
which are obtained in a frequentist approach. We find 
that even if = 0 is excluded at high significance in 
a Bayesian approach this is no longer the case from the 
frequentist point of view. In other words the probability 
that a model with vanishing Ha leads to the present-day 
observed CMB and LSS data is not exceedingly small. 
We then study how the results are modihed if we allow 
for general isocurvature contributions to the initial condi¬ 
tions SSIl- In this first study of the matter power 
spectrum from general isocurvature modes we discover 
that a COBE-normalized matter power spectrum repro¬ 
duces the observed amplitude only if it is highly domi¬ 
nated by the adiabatic component. Hence the isocurva¬ 
ture modes cannot contribute significantly to the matter 
power spectrum and do not lead to a degeneracy in the 
initial conditions for the matter power spectrum when 
combined with CMB data. This is the main result of our 
paper. 

The paper is organized as follows: In the next section 
we discuss the setup for our analysis, the space of cosmic 
parameters and of initial conditions, and we recall the dif¬ 
ference between Bayesian and frequentist approach. In 
Sec. 3 we present the results for adiabatic and for mixed 
(adiabatic and isocurvature) perturbations. Sec. 4 is ded¬ 
icated to the conclusions. 


II. ANALYSIS SETUP 
A. Cosmological parameters 

As it has been discussed in the literature, the recent 
data-sets, BOOMERanG §, MAXIMA |T|, DASI ||, 
VSA 1^, CBI and Archeops Q are in very good 
agreement up to the third peak in the angular power 
spectrum of CMB anisotropies, £ ^ 1000. In our anal¬ 
ysis we therefore use the COBE data [|^ (7 points ex¬ 
cluding the quadrupole) for the £ region 3 < £ < 20 
and the BOOMERanG data to cover the higher £ part of 
the spectrum (19 points in the range 100 < £ < 1000). 
Since Archeops has the smallest error bars in the region 
of the first acoustic peak, we also include this data-set 
(16 points in the range 15 < ^ < 350). Including any of 
the other mentioned data does not influence our results 
signihcantly. The BOOMERanG and Archeops absolute 
calibration errors (10% and 7% at Icr, respectively) as 
well as the uncertainty of the BOOMERanG beam size 


are included as additional Gaussian parameters, and are 
maximized over. We make use of the Archeops window 
functions found in Ref. |Q, while for BOOMERanG a 
top-hat window is assumed. 

For the matter power spectrum, we use the galaxy- 
galaxy power spectrum from the 2dF data which is ob¬ 
tained from the redshift of about 10^ galaxies [^. We in¬ 
clude only the 22 decorrelated points in the linear regime, 
i.e., in the range 0.017 < k < 0.314 [h Mpc“^], and the 
window functions of Ref. [pot which can be found on the 
homepage of M. Tegmark |39f|. 

Our grid of models is restricted to flat universes and we 
assume purely scalar perturbations. Since the goal of this 
paper is more to make a conceptual point than to con¬ 
sider the most generic model, we fix the baryon density 
to the BBN preferred value = Wb = 0.020 |2^. We 
investigate the following 3 dimensional grid in the space 
of cosmological parameters: 0.80 < ns (0.05) < 1.20, 
0.35 < h (0.025) < 1.00, 0.00 < Ha (0.05) < 0.95, where 
ns is the scalar spectral index and the numbers in paren¬ 
thesis give the step size we use. The total matter content 
Hm = He -l- Hb is Hni = 1 — Ha, and He indicates the 
cold dark matter contribution. For all models the optical 
depth of reionization is r = 0 and we have 3 families of 
massless neutrinos. For each point in the grid we com¬ 
pute the ten GMB and matter power spectra, one for 
each independent set of initial conditions (see Sec. HB 
below). 


B. Allowing for isocurvature modes 

We enlarge the space of models by including all pos¬ 
sible isocurvature modes. As it has been argued in 
Ref. ||l^, generic initial conditions for a fluid consisting 
of photons, neutrinos, baryons and dark matter allow 
for five relevant modes, i.e., modes which remain regu¬ 
lar when going backwards in time. These are the usual 
adiabatic mode (AD), the cold dark matter isocurvature 
mode (GI), the baryon isocurvature mode (BI), the neu¬ 
trino isocurvature density (NID) and neutrino isocurva¬ 
ture velocity (NIV) modes. The CMB and matter power 
spectra from the cold dark matter and the baryon isocur¬ 
vature modes are identical (see the argument given in 
Ref. [^) and therefore from now on we will just con¬ 
sider one of those, namely the Cl mode. We assume that 
all these four modes are present with arbitrary initial 
amplitude and arbitrary correlation or anti-correlation. 
The only requirement is that their superposition must 
be a positive quantity, since the C(S and matter power 
spectrum are quadratic and thus positive observables. 
For simplicity we restrict ourselves to the case where all 
modes have the same spectral index. Initial conditions 
are then described by the spectral index ns and a posi¬ 
tive semi-definite 4x4 matrix, which amounts to eleven 
parameters instead of two in the case of pure AD initial 
conditions. More details can be found in Ref. ||l^. For 
the search among initial conditions we use the simplex 
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method described in Ref. |T^, with the following modi¬ 
fication. We find it convenient to express the matrix A 
describing the initial conditions as 

A = UDU^, (1) 

where A G Vn, U G *SO„, D = diag((ii, ^ 2 , ■ • ■, dn) and 
di > 0, i G {l,2,...,n}. Here Vn denotes the space of 
nxn real, positive semi-definite, symmetric matrices and 
SOn is the space of n x n real, orthogonal matrices with 
det = 1. As explained above, here we take n = 4. We 
can write U as an exponentiated linear combination of 
generators of SOn- 


U = exp 


' {rd-n)l1 




( 2 ) 


with 


Hi 


/ 0 1 0 ...\ 
-1 0 0 ... 

0 0 0 ... 

V ; ; ; ■.) 


(3) 


and so on, with —7r/2 < Ui < 7r/2, i G {1,2,..., (n^ — 
n)/2}. In analogy to the Euler angles in three dimen¬ 
sions, we can re-parameterize U in the form 


{n'^-n)l2 

U= n exp(i/>,H,), (4) 

i=l 


with some other coefficients —tt/ 2 < tpi < ttI 2, i G 
{1,2,..., (n^ — n)/2}, whose functional relation with the 
ads does not matter. The diagonal matrix D can be 
written as 


D = diag(tan(6»i),... ,tan(0„)), (5) 

with 0 < < 7r/2, for i G {l,2,...,n}. In this 

way, the space of initial conditions for n modes is ef¬ 
ficiently parameterized by the (n^ -|- n)/2 angles 
In our case, n = 4 and the initial conditions are de¬ 
scribed by the ten dimensional hypercube in the variables 
(01,..., 04 , i/ji,..., i/jg)- This is of particular importance 
for the numerical search in the parameter space. One can 
then go back to the explicit form of A using Eqs. 
and (1). 

For a given initial condition determined by a pos¬ 
itive semi-definite matrix A and a spectral index ns 
we quantify the isocurvature contribution to the CMB 
anisotropies by the parameter /3 defined as 

y: (Wf+iHcy"), 

^ X = CI,NIV,NID 

y = AD,CI,NIV,NID 

where the average <, > is taken in the £ range of interest, 
in our case 3 < £ < 1000, and where stands for 

the auto-correlator of the CMB anisotropies with initial 
conditions X. 


C. Bayesian or frequentist? 


For the sake of clarity, we briefly recall two possible 
points of view which one can take when doing data analy¬ 
sis, the Bayesian and the frequentist approach, and high¬ 
light their difference. More details can be found, e.g., in 
Refs. 1^, 1^, Another possibility is based on Markov 
Chain Monte Carlo techniques, which we do not discuss 
here; see instead and references therein. 

When fitting experimental data, we minimize a over 
the parameters which we are not interested in. This 
procedure is equivalent to marginalization if the random 
variables are Gaussian distributed. The Maximum Like¬ 
lihood (ML) principle states that the best estimate for 
the unknown parameters 0 is the one which maximizes 
the likelihood function: 


C = Coexp{-x^/2). (7) 


We then draw Icr, 2cr and 3cr likelihood contours around 
the ML point, i.e., the one for which the x^ is minimal in 
our grid of models. The likelihood contours are defined to 
be A = — Xml = 2.30, 6.18, 11.83 away from the ML 

value for the joint likelihood in two parameters, A = 1, 
4, 9 for the likelihood in only one parameter. This is the 
Bayesian approach: in a somewhat fuzzy way, likelihood 
intervals measure our degree of belief that the particular 
set of observations used in the analysis is generated by a 
parameter set belonging to the specified interval. In this 
case, one implicitly accepts the ML point in parameter 
space as the true value, while points which are further 
away from it are less and less “likely” to have generated 
the measurements. This is the content of Bayes’ The¬ 
orem, which allows us to interpret the joint conditional 
probability £(x|0) of measuring x for a fixed set of pa¬ 
rameters 0 as the inverse probability P(0|x) for the 
value of 0 given the measurements x. 

On the other hand, in the frequentist approach one 
asks a different question: What is the probability of ob¬ 
taining the experimental data at hand, if the Universe 
has some given cosmological parameters, e.g., a vanish¬ 
ing cosmological constant? Clearly, if we want to answer 
the question whether a certain set of experimental data 
forces a non-vanishing cosmological constant, this is ac¬ 
tually the correct question to ask. To the extent to which 
the CiS can be approximated as Gaussian variables, the 
quantity x^ is distributed according to a chi-square prob¬ 
ability distribution with F = N — M degrees of free¬ 
dom (dof), which we denote by where N is the 

number of independent (uncorrelated) experimental data 
points and M is the number of htted parameters. From 
the distribution one can readily estimate confidence 
intervals. For a given parameter set 0 with chi-square 
X^ the probability that the observed chi-square will be 
larger than the actual value by chance fluctuations is 



P^^\x)dx = 1 — 7 . 


( 8 ) 
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In other words, if the measurements could be repeated 
many times in different realizations of our universe, the 
estimated confidence interval would asymptotically in¬ 
clude the true value of the parameters 1007 % of fho time. 

It is customary in cosmological parameter estimation 
to present likelihood plots drawn using the Bayesian ap¬ 
proach. It is misleading that such Bayesian contours are 
usually called “confidence contours”, which properly des¬ 
ignate frequentist contours. Likelihood (Bayesian) con¬ 
tours are usually much tighter than the confidence con¬ 
tours drawn from the frequentist point of view. This 
is a consequence of the ML point having often a 
much smaller than 1, because the data-sets are highly 
consistent with each other and also because usually not 
all points are completely independent. If we consider the 
usual situation in which likelihood contours are drawn in 
a two dimensional plane with all other parameters maxi¬ 
mized, the frequentist approach is more conservative than 
the Bayesian one. This is because in general, for reason¬ 
ably good ML values x^i^<0{F) and F > 2, 

nCO 

P^^\x)dx = / P^^\x)dx (9) 

hi 

only for When looking at likelihood contours 

one should thus keep in mind that a point more than 
say 3cr away from the ML point is not necessarily ruled 
out by data, as we shall show below. In order to es¬ 
tablish this, one has to look at confidence contours, i.e., 
ask the frequentist’s question. In the following, the term 
“likelihood contours” will refer to contours drawn in the 
Bayesian approach, while the term “confidence contours” 
will be reserved for contours coming form the frequentist 
point of view. 

III. RESULTS 
A. Adiabatic perturbations 

We first fit CMB data only [N = 42) by maximizing 
M = 7 parameters, i.e., the BOOMERanG and Archeops 
calibration errors, BOOMERanG beam size error, ng, h, 
Oa and the overall amplitude of the adiabatic spectrum, 
and we find (Bayesian likelihood intervals on Oa alone): 

= 0.801)J;3^ at 2cr and to;8o at 3cr. (10) 

The asymmetry in the intervals arises because the value 
of Oa for our ML model is relatively large. One could 
achieve a better precision in determining the ML value 
of Oa by using a finer grid and varying Wb as well, which 
has extensively been done in the literature and is not the 
scope of this work. Moreover, the position of the acous¬ 
tic peaks in GMB anisotropies is mainly sensitive to the 
age of the universe at recombination, which depends only 
on and to the angular diameter distance, which 

depends on Om, Oa and the curvature of the universe. 
When the universe is flat, the angular diameter distance 


is weakly dependent on the relative amounts of flm and 
ft A as soon as Ha is not too large (see e.g. Ref. |^ ). 
Hence, one can achieve a sufficiently low value of 
either via a large cosmological constant or via a very low 
Hubble parameter, /i<0.45. 

We now include the matter power spectrum Pm, as¬ 
suming Pm = 6^Pg, where Pg is the observed galaxy 
power spectrum and b some unknown bias factor (as¬ 
sumed to be scale independent), over which we maxi¬ 
mize. Inclusion of this data in the analysis breaks the 
Ha, h degeneracy, since Pm is mainly sensitive to the 
shape parameter T = Hm/i. We therefore obtain signifi¬ 
cantly tighter overall likelihood intervals for Ha: 

Ha = at 2cr and at 3a. (11) 

We plot joint likelihood contours for Ha, h with AD 
initial conditions in Fig. From the Bayesian analy- 



FIG. 1: Joint likelihood contours (Bayesian), with CMB only 
(solid lines, showing la, 2a, 3a contours) and CMB-I-LSS 
(filled) for purely adiabatic initial conditions. 

sis, one concludes that GMB and LSS together require a 
non-zero cosmological constant at very high significance, 
more than 7a for the points in our grid! Note that the 
ML point has a reduced chi-square Xf =56 = 0.59, signif¬ 
icantly less than 1. 

The frequentist analysis, however, excludes a much 
smaller region of parameter space (Fig. |^). The frequen¬ 
tist contours must be drawn for the effective number of 
dof, i.e., using the number of effectively independent data 
points. We can therefore roughly take into account a 10% 
correlation, which is the maximum correlation between 
data points given in 0 ,@, by replacing F by the effec¬ 
tive number of dof, Pes = 0.9N — M, and rounding to 
the next larger integer (to be conservative). One could 
argue that the BOOMERanG and Archeops data points 
are not completely independent, since BOOMERanG ob¬ 
served a portion of the same sky patch as measured by 
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Archeops. This possible correlation is difficult to quan¬ 
tify, but should not be too important since the sky por¬ 
tion observed by Archeops is a factor of 10 larger than 
BOOMERanG’s and therefore we ignore it here. Fig. ^ 
is drawn with Fgg = 31 for CMB alone and Fgs = 50 for 
CMB-I-LSS, but we have checked that our results do not 
change much if we use a 5% correlation. It is interesting 
to note that there are regions in Fig. which are excluded 
with a certain confidence by CMB data alone but are no 
longer excluded at the same confidence when we include 
LSS data. In other words, it would seem that taking into 
account more data and therefore more knowledge about 
the universe, does not systematically exclude more mod¬ 
els, i.e., the CMB-I-LSS contours are not always contained 
in the CMB alone contours. This apparent contradiction 
vanishes when one realizes that the conhdence limits on, 
e.g., Ha alone in the frequentist approach are just the 
projection of the confidence contours of Fig. ||on the Qa 
axis. One can readily verify in Fig. ^ that the confidence 
limits for the combined data-set are always smaller than 
the ones for CMB data alone. There are points with 
Oa = 0 and h ~ 0.40 which are still compatible within 
2(7 with both LSS and CMB data, at the price of push¬ 
ing somewhat the other parameters. In the best fit with 
Oa = 0 shown in Fig. ||, one has to live with a red spec¬ 
tral index ns = 0.80. Furthermore, the calibration of the 
BOOMERanG and Archeops data points is reduced in 
this fit by 34% and 26%, respectively, i.e., more than 3 
times the quoted la systematic error. 



FIG. 2: Confidence contours (frequentist) with CMB only 
(solid lines, la, 2a, 3a contours and FeS = 31) and 
CMB-I-LSS (filled, Fefi = 50) for purely adiabatic initial con¬ 
ditions. 

In both cases, it is clear that one can exploit the SIa, h 
degeneracy to fit CMB data alone with a model having 
Ha = 0. For a flat universe like the one we are con¬ 
sidering, one has then to go to a much smaller value of 
the Hubble parameter than the one indicated by other 




FIG. 3: Best fit with Qa ~ 0 and purely AD initial condi¬ 
tions, compatible with CMB and LSS data within 2a con¬ 
fidence level. In the lower panel, only the 2dF data points 
left of the vertical, dotted line — i.e., in the linear region — 
have been included in the analysis. Note the low first acoustic 
peak due to the joint effect of the red spectral index and of 
the absence of early ISW effect. In this fit, the calibration 
of BOOMERanG (red errorbars) and Archeops (blue error- 
bars) has been reduced by 34% and 26%, respectively. This is 
more than 3 times the quoted la calibration errors for both 
experiments. To appreciate the difference, we plot the non 
recalibrated value of the BOOMERanG and Archeops data 
points as light blue and magenta crosses, respectively. In the 
upper panel, green errorbars are the COBE measurements. 
Even though the fit is “by eye” very good, it seems highly 
unlikely that the calibration error is so large. 


measurements, most notably the HST Key Project 
which gives h = 0.72 ± 0.08. The LSS data are mainly 
sensitive to the shape parameter T ~ 0.2. Hence LSS 
with = 1.0 would require an even lower value of h 
which is not compatible with CMB. Therefore inclusion 
of LSS data tends to exclude any flat model without a 
cosmological constant. Summing up, for purely adiabatic 
initial conditions the Bayesian approach gives very strong 
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support to r^A 7 ^ 0; in the more conservative frequentist 
point of view, while Qa 7^ 0 cannot be excluded with 
very high confidence, present LSS and CMB data start to 
be incompatible with a flat universe with vanishing cos¬ 
mological constant. These conclusions are in qualitative 
agreement with previous works |3^, |2[ jsj, HI . 

In the next section we investigate the stability of those 
well known results with respect to inclusion of non- 
adiabatic initial conditions. 


B. Isocurvature modes 

We now enlarge the space of models by including all 
possible isocurvature modes with arbitrary correlations 
among themselves and the adiabatic mode as described 
in the previous section. We first consider CMB data only 
and maximize over initial conditions. The number of pa¬ 
rameters increases by nine and the number of dof de¬ 
creases correspondingly with respect to the purely AD 
case considered above. Likelihood (Bayesian, see Fig. 
and confidence (frequentist, see Fig. ^ contours widen up 
somewhat along the degeneracy line. The enlargement is 
less dramatic than for other parameter choices, see, e.g.. 
Ref. |l^ where the degeneracy in ujh,h was analyzed. 
This is partially due to our prior of flatness which reduces 
the space of models to the ones which are almost degener¬ 
ate in the angular diameter distance. Most of our models 
have the first acoustic peak of the adiabatic mode already 
in the region preferred by experiments. Hence in most of 
the fits isocurvature modes play a modest role, especially 
in the parameter regions with large Da, h (cf Fig. ^ and 
the discussion below). Nevertheless, because of the Da, h 
degeneracy, even a modest widening of the contours along 
the degeneracy line results in an important worsening of 
the likelihood limits. The ML point does not depart very 
much from the purely adiabatic case, but now we cannot 
constrain Da at more than Icr (Bayesian, CMB only): 


Da = at la, (12) 

and no limits for 0.0 < Da < 0.95 at higher confidence. 

In Fig. II we plot the dark matter power spectra of the 
different auto- (upper panel) and cross-correlators (lower 
panel) for a concordance model. The norm of each pure 
mode (AD, Cl, NID, NIV) is chosen such that the cor¬ 
responding CMB power spectrum is COBE-normalized. 
The cross-correlators are normalized according to totally 
correlated spectra, i.e. 

^(x.Y) = V ^xAy/2 , (13) 

where A(x,y) denotes the norm of the cross-correlator 
between the modes X,V and Ax the norm of the pure 
mode X. The CMB power spectrum for this set of cos¬ 
mological parameters can be found in Ref. . A crucial 
result is that the COBE-normalized amplitude of the AD 
matter power spectrum is nearly two orders of magnitude 



FIG. 4: Joint likelihood contours (Bayesian) with general 
isocurvature initial conditions, with CMB only (solid lines, 
Ict, 2ct, 3a) and CMB-l-LSS (filled). The disconnected Icr 
region is an artificial feature due to the grid resolution. 


larger than the isocurvature contribution. The main rea¬ 
son for this is the amplitude of the Sachs Wolfe plateau 
which is about for adiabatic perturbations and 2'!' for 
isocurvature perturbations. Here T' is the gravitational 
potential. This difference of a factor of about 36 in the 
power spectrum on large scales is clearly visible in the 
comparison of Pad and Pci (the difference increases at 
smaller scales). The case of the neutrino modes is even 
worse since they start up with vanishing dark matter per¬ 
turbations. That the CDM isocurvature matter power 
spectrum is much lower than the adiabatic one has been 
known for some time (see e.g. Ref. js^]). However, it was 
not recognized before that the same holds true for the 
neutrino isocurvature matter power spectra as well, and 
- more importantly - that this leads to a way to break 
the strong degeneracy among initial conditions which is 
present in the CMB power spectrum alone. 

In an analysis with general initial conditions includ¬ 
ing LSS data only we obtain very broad likelihood and 
confidence contours which exclude only the lower right 
corner of the (Da,/i) plane. In contrast to the CMB 
power spectrum, the matter power spectrum can be fit¬ 
ted with extremely high adiabatic and isocurvature con¬ 
tributions, which are then typically cancelled by large 
anti-correlations between the spectra. This behavior is 
exemplified for a model with general IC and Da = 0.70, 
h = 0.65, ns = 1.0 in Fig. ^ The best fits with LSS 
data only are dominated by large isocurvature cross¬ 
correlations. Clearly, the resulting CMB power spectrum 
is highly inconsistent with the COBE data. Hence such 
“bizarre” possibilities are immediately ruled out once we 
include CMB data. Conversely, moderate isocurvature 
contributions can help fitting the CMB data, and do not 
influence the matter power spectrum, which is completely 
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FIG. 5: Dark matter power spectra of the different auto- 
(upper panel) and cross-correlators (lower panel) for a con¬ 
cordance model with Qa ~ 0.70, h = 0.65, ns = 1.0, 
u!b = 0.020, with the corresponding CMB power spectrum 
COBE-normalized (see the text for details). The color 
codes are as follows: in the upper panel, AD: solid/black 
line. Cl: dotted/green line, NID: short-dashed/red line, NIV: 
long-dashed/blue line; in the lower panel, AD: solid/black 
line (for comparison), < AD, Cl >: long-dashed/magenta 
line, < AD, NID >: dotted/green line, < AD, NIV >: 
short-dashed/red line, < Cl, NID >: dot-short dashed/blue 
line, < Cl, NIV >: dot-long dashed/light-blue line, and 
< NID, NIV >: solid/yellow line. The adiabatic mode is by 
far dominant over all others. 

dominated by the AD mode alone. Combining CMB and 
LSS data (see Fig. ^ we find now (Bayesian, mixed IC): 

= 0.65l^;25 at 2cr and at 3cr. (14) 



FIG. 6: Concordance model fit with general IC and LSS 
data only. The total spectrum (solid/black) is the re¬ 
sult of a large cancellation of the purely AD part (long- 
dashed/red) by the large, negative sum of the various cor¬ 
relators (dotted/magenta, plotted in absolute value). The 
short-dashed/green curve is the sum of the three pure isocur¬ 
vature modes. Cl, NID and NIV. Note that the resulting total 
spectrum is less than one tenth of the purely adiabatic part. 

combination of CMB and LSS measurements turn out to 
be an ideal tool to constrain the isocurvature contribu¬ 
tion to the initial conditions. 



FIG. 7: Confidence contours (frequentist) with general isocur¬ 
vature initial conditions with CMB only (solid lines, Icr, 2a, 
Sct contours Fes = 22) and CMB-I-LSS (filled, Fbb = 41). 


The likelihood limits are larger than for the purely adia¬ 
batic case but it is interesting that the Bayesian analysis 
still excludes Ha = 0 at more than 3cr even with general 
initial conditions, for the class of models considered here. 
Because of the above explained reason, the widening of 
the limits is not as drastic as one might fear. Therefore, 


From the frequentist point of view, one notices that the 
region in the Da, h plane which is incompatible with data 
at more than 3tT is nearly independent on the choice of 
initial conditions (compare Fig. ^and Fig. Enlarging 
the space of initial conditions seemingly does not have a 
relevant benefit on fitting present-day data with or with- 




































out a cosmological constant. The reason for this is that 
the (COBE-normalized) matter power spectrum is dom¬ 
inated by its adiabatic component and therefore the re¬ 
quirement ~ 0.2 remains valid. In Fig. ^ we plot the 
best fit model with general initial conditions and IIa = 0. 
We summarize our likelihood and confidence intervals on 
IIa (this parameter only) in Table 1. 




FIG. 8: Best fit with general IC and Q.a = 0. As for purely 
AD, even with general IC the absence of the cosmological 
constant suppresses in an important way the height of the 
first peak. In both panels we plot the best total spectrum 
(solid/black), the purely AD contribution (long-dashed/red), 
the sum of the pure isocurvature modes (short-dashed/green) 
and the sum of the correlators (dotted/magenta, multiplied 
by —1 in the upper panel and in absolute value in the lower 
panel). The matter power spectrum is completely dominated 
by the AD mode, while the correlators play an important 
role in cancelling unwanted contributions in the CMB power 
spectrum at the level of the first peak and especially in the 
COBE region. For this model we have /3 = 0.39, while the 
BOOMERanG and Archeops calibrations are reduced by 28% 
and 12%, respectively. 

In Fig. 1^ we plot the isocurvature contribution to the 
best fit models with CMB and LSS in terms of the pa¬ 


rameter f3 defined in Eq. (||). The best fit with IIa = 0 
has an isocurvature contribution of about 40%. 



FIG. 9: Isocurvature content 0.0 < /3 < 1.0 of best fit 
models with CMB and LSS data. The contours are for 
(3 = 0.20, 0.40, 0.60, 0.80 from the center to the outside. 

It would be of great interest to investigate whether the 
result FIa 7 ^ 0 is robust with respect to addition of general 
initial conditions in an open universe 

IV. CONCLUSIONS 

The conclusions of this work are threefold. The first 
one is not new, but seems to be dangerously forgotten 
in recent cosmological parameters estimation literature: 
namely that likelihood contours cannot be used as “ex¬ 
clusion plots”. The latter are usually substantially wider, 
less stringent. A more rigorous possibility are frequentist 
probabilities, which however suffer from the dependence 
on the number of really independent measurements which 
is often very difficult to come by. 

Secondly, we have found that in COBE-normalized 
fluctuations, the matter power spectrum has negligible 
isocurvature contributions and is essentially given by the 
adiabatic mode. Hence the shape of the observed matter 
power spectrum still requires Om/i ~ 0.2, independent of 
the choice of initial conditions. Due to this behavior, the 
condition F2 = OA+Om = 1 requires either a cosmological 
constant or a very small value for the Hubble parameter, 
independently from the isocurvature contribution to the 
initial conditions. 

The third conclusion from our work are the following 
results for the presence of a cosmological constant: For 
flat models, a likelihood (Bayesian) analysis strongly fa¬ 
vors a non-vanishing cosmological constant. Even if we 
allow for isocurvature contributions with arbitrary corre¬ 
lations, a vanishing cosmological constant is still excluded 
at more than 3 (t. If we would allow for open models, a 
significant contribution from the NIV mode which has 
the first acoustic peak at £ = 170 in flat models, possibly 
could at the same time give a good fit to CMB data and 
allow for the observed shape parameter T with a reason¬ 
able value of h. For technical reasons we shall study this 
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TABLE I: Results for the likelihood (Bayesian) and confidence (frequentist) intervals for alone (all other parameters 
maximized). A bar, —, indicates that at the given likelihood/confidence level the analysis cannot constraint Ha in the range 
0.0 < Ha < 0.95. Where the quoted interval is smaller than our grid resolution, an interpolation between models has been 
used. 


AD only 



Bayesian 


Frequentist 




Data-set 

SIa 

la 

2a 3a 


la 

2a 

3f7 

F 

xVf 

CMB“-|-flatness 

0.80 

+ 0.08 

- 0.08 

+ 0.10 + 0.12 
- 0.35 


< 0.93 

- 

- 

35 

0.58 

CMB“ -hLSS‘-hfiatness 

0.70 

+ 0.05 

- 0.05 

+ 0.13 + 0.15 
- 0.17 - 0.27 

0.15 < Da < 0.90 

< 0.92 

< 0.92 

56 

0.59 

General IC 

CMB“-|-flatness 

0.85 

+ 0.05 

- 0.35 

- 


- 

- 

- 

26 

0.74 

CMB“ -hLSS''-hfiatness 

0.65 

+ 0.15 

- 0.10 

+ 0.22 + 0.25 
- 0.25 - 0.48 


< 0.90 

< 0.92 

< 0.95 

47 

0.67 


“ COBE, BOOMERanG and Archeops data. 

^ 2dF data. 

Likelihood interval. 

Region not excluded by data with given confidence. 


case in a forthcoming paper . 

The situation changes considerably in the frequen¬ 
tist approach. There, even for purely adiabatic models, 
Ha = 0 is still within 3cr for a value of h < 0.48 which is 
marginally defendable. The conclusion does not change 
very much when we allow for generic initial conditions. 
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